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Qj [ cone. In this approach, the large-distance dynamics is incorporated in a set of pion wave 

' functions. We take into account two-particle and three-particle wave functions of twist 

2, 3 and 4. The resulting values of the coupling constants are qd^d-k = 12.5 ± 1 and 
; gB*Bn = 29 ± 3. From this we predict the partial width T{D*+ D^tt+) = 32 ± 5 keV . 

■ We also discuss the soft-pion limit of the sum rules which is equivalent to the external 

axial field approach employed in earlier calculations. Furthermore, using gs* b-k and 
gD*Dn the pole dominance model for the B ^ tt and D ^ n semileptonic form factors 
is compared with the direct calculation of these form factors in the same framework of 
light-cone sum rules. 
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1 Introduction 



The extraction of fundamental parameters from data on heavy flavoured hadrons inevitably 
requires some information about the physics at large distances. Numerous theoretical studies 
have been devoted to making this extraction as reliable as possible. While the inclusive B 
and D decays appear to be the most clean reactions theoretically, exclusive decays are often 
much easier to measure experimentally. However, for their interpretation one needs accurate 
estimates of decay form factors and other hadronic matrix elements. In the exceptional case 
B Deu, the form factor at zero recoil can be calculated in the heavy quark limit |^. 
In most other important cases, one has to rely on less rigorous nonperturbative approaches. 
Among those, QCD sum rules 0] have proved to be particularly powerful. 

In this paper we employ sum rule methods in order to calculate the 0*071 and B*Bt^ 
coupling constants. These couplings are interesting for several reasons. In particular, they 
determine the normalization of the heavy-to-light form factors D ^ tt and B ^ n near zero 
pionic recoil, where the D*- and B*- poles are believed to dominate. For further discussion 
one may consult refs. 0, ^. Recently, it has been argued that in the combined heavy quark 
and chiral limit vector meson dominance becomes even exact 0. As noted in refs. J^, §[, B* 
dominance is also compatible with the dependence of the B tt form factor on the momentum 
transfer p predicted by QCD sum rules at low values of p'^. A similar conclusion is reached in 
ref. concerning D* dominance in the D ^ tt form factor. In addition, the D*Dti coupling 
can be directly measured in the decay D* Dtt and thus provides one more independent test 
of the sum rule approach. 

Calculations of couplings of heavy mesons to a pion have already been undertaken several 
times in the framework of QCD sum rules. Unfortunately, the sum rules obtained in refs. 
[10-13] differ in nonleading terms and, to some extent, also in numerical results. Here, we 
suggest an alternative method known as QCD sum rules on the light-cone. In this approach, 
the ideas of duality and matching between parton and hadron descriptions, intrinsic to the 
QCD sum rules, are combined with the specific operator product expansion (OPE) techniques 



used to study hard exclusive processes in QCD [14, 15|. In contrast to the conventional sum 
rules based on the Wilson OPE of the T-product of currents at small distances, one considers 
expansions near the light-cone in terms of nonlocal operators, the matrix elements of which 
define hadron wave functions of increasing twist. As one advantage, this formulation allows to 
incorporate additional information about the Euclidean asymptotics of correlation functions 
in QCD for arbitrary external momenta. These new features are related to the (approximate) 
conformal invariance of QCD and are coded in the hadron wave functions. Many of the 
theoretical results obtained in the context of exclusive processes (see e.g. ref. [|l6l) are very 
useful in the present context as well. In turn, we will see that heavy-flavour decays can provide 
valuable constraints on the wave functions. 

Previous applications of light-cone sum rules include calculations of the amplitude of the 
radiative decay S ^ p7 [O, the nucleon magnetic moments [0, the strong couplings g^NN 



and (7ptj7r ^rm factors of semileptonic and radiative B- and D-meson decays [8,19-21], the 
pion form factor at intermediate momentum transfers [^], and the ttAj* form factors [^]. In 
all these cases the results are encouraging. 

The light-cone sum rule for the coupling of heavy mesons to a pion is the principal result 
of the present paper which is organized as follows. In Sect. 2 we discuss possible strategies 
in constructing sum rules for coupling constants and explain the concept of light-cone sum 
rules. The derivation of the sum rule for the B*Bit and D*Dtt couplings is then completed 



1 



in Sect. 3, taking into account the pion two- and three-particle wave functions up to twist 4. 
Sect. 4 is devoted to a detailed numerical analysis. In Sect. 5 we show that in a simplified case, 
putting the external momenta in the correlation function equal to each other and performing 
a Borel transformation in one momentum instead of two, we obtain the sum rule proposed 
previously in refs. [10-13]. We demonstrate that despite the slightly different terminology of 
these papers the sum rules must coincide with each other, and elaborate on possible subtleties 
in these earlier calculations. Furthermore, in Sect. 6 using our results on the B*Btt and D*Dtt 
coupling constants, we confront the pole model for the heavy-to-light form factors B ^ tt and 
D —>■ TT with a direct calculation of these form factors in the same framework of light-cone sum 
rules following ref. [§. A comprehensive comparison of our results on the 0*071 and B*Btt 
couplings with other estimates and our conclusions are presented in Sect. 7. 

Technical details are collected in two Appendices. Appendix A summarizes the relevant 
features of the pion wave functions and specifies the input in our numerical calculations. In 
Appendix B we derive a simple rule how to subtract the contribution from excited resonances 
and continuum states in the sum rule. 



2 Light-cone versus short-distance expansion 

For definiteness, we focus on the D*D7i coupling defined by the on-mass-shell matrix element 

{D*^{p)n-{q)\D%p + q)) = -gD*D.q,e\ (1) 

where the momentum assignment is specified in brackets and is the polarization vector of 
the D* . The couplings for the different charge states are related by isospin symmetry: 

gD*D-K = 9d*+D07t+ = -V2gD* + D+nO = V2gD*ODOnO = -QD'OD+t,- ■ (2) 

Most of what is said below applies equally to B*B7i couplings. The corresponding relations 
are obtained by the obvious replacements c b, D* — * B* and D ^ B. 

Following the general strategy of QCD sum rules, we want to obtain quantitative esti- 
mates for qd* /jtt by matching the representations of a suitable correlation function in terms of 
hadronic and quark-gluon degrees of freedom. For this purpose, we choose 

F^{p,q) = I I rfWP^(7r-(g)|T{J(x)7^c(x),c(0)z75w(0)}|0) . (3) 



To the best of our knowledge, the study of correlation functions with the T-product of currents 
sandwiched between the vacuum and one-pion state was first suggested in ref. . 

With the pion being on mass-shell, = m^, the correlation function (^) depends on two 
invariants, p"^ and (p + q)"^- Throughout the paper we set 171^^ = 0. The contribution of interest 
is the one having poles in and {p + q^ '■ 

TP ( \ rnlniD-fofD-gD-Dn 1 mjy 



m 



D* 



Obviously, this term stems from the ground states in the {dc) and {cu) channels. To derive 
eq. we have made use of eq. (P and the decay constants fn and fo* defined by the matrix 
elements 

(D|d75w|0) = (5) 
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and 



(0 I d-i^c I D*) = mD*fD*e^, , (6) 

respectively. 

The main theoretical task is the calculation of the correlation function (|^) in QCD. This 
problem can be solved in the Euclidean region where both virtualities and {p + qY are 
negative and large, so that the charm quark is sufficiently far off-shell. Substituting, as a first 
approximation, the free c-quark propagator 

(0|T{cWc(0)}|0) = ,£fi(x) = / ^e-'-ijtBL (7) 



into eq. (H) one readily obtains 



(27r)4(m2 - P) 
+ r(7r(g)|J(x)7^7.75«(0)|0) 

Diagramatically, this contribution is depicted in Fig. la. Applying the short- distance expan- 
sion (SDE) in terms of local operators to the first matrix element of eq. (^, 

d{x)-f^l,>u{^) = E -^m(D ■xri.l.ui^) , (9) 

one has after integration over x and k : 

F.iP^l) = ^^^^flr^P^uMnq,. (10) 
m-^ — p'^ [mi — p^)"- 

where 

{7r{q)\d Da,Da2 ■■■ Da„ 7m75m|0) = {iTqt,qa^qa2-qa„Mn + ••• , 

D being the covariant derivative, has been used. One immediately encounters the following 
problem. If the ratio 

e = 2{p ■ q)/{ml - = ((p + qf - p2)/(^2 _ ^^^^ 

is finite one must keep an infinite series of local operators in eq. (]T0|). All these operators give 
contributions of the same order in the heavy quark propagator l/(m^ — p^), differing only by 
powers of the dimensionless parameter ^ 0. Therefore, SDE of eq. (|^) is useful only if ^ ^0, 
i.e. if p2 ~ (p -|- qY or, equivalently, g ~ 0. Under this condition, the series in eq. (|T0|) can 
be truncated after a few terms involving only a small number of unknown matrix elements 
Mn- However, for general momenta with p^ ^ {p + qY oiie has to sum up the infinite series of 
matrix elements of local operators in some way. 

This formidable task is solved by using the techniques developed for hard exclusive pro- 



cesses in QCD |1T5| , |T6[. We illustrate the solution for the correlation function 

I j d'xe'^%n\q)\T{q{x)^,Qq{x)A{^)^Mqm\^) = e,uc.pp''q^F\p\{p + q)^ , (12) 

-•^This feature is also observed in deep inelastic scattering, with the variables {Q^, v, x} playing the role of 
{— p^, p ■ q, ^} . As well known, there one applies an expansion near the light-cone in terms of operators of 
increasing twist, rather than of increasing dimension. 



3 



which is similar to eq. and defines the form factor of the couphng of a pion to a pair of 
virtual photons |]T6|, In eq. (p!2D, Q is a matrix of electromagnetic charges, and g is a row 



vector composed of the up and down quark fiavours. As well-known |Tj], for sufficiently virtual 
photons, — oo and {p + g)^ — oo, this form factor can be calculated in perturbative 

QCD. The principal result reads 

F-(/.(p + #) = f„/'i?^, f„ = i^. (13) 



(p + uqY 

with calculable radiative corrections, and with power corrections suppressed by the photon 
virtualities. Here, <fTr{u) is the pion wave function of leading twist, defined by the following 
matrix element of a nonlocal operator on the light-cone = : 

(7r(g)|J(x)7^75^(0)|0) = -iq^U f' due'^^'^^Mu) . (14) 

Jo 

Physically, (p-,^ represents the distribution in the fraction of the light-cone momentum go + Qs 
of the pion carried by a constituent quark. Note the normalization of v^tt to unity following 
from eq. (|l^) for x = . 

Let us first concentrate on the form factor (0) at (almost) equal photon virtualities, i.e. 
at ^ = (2p • q)/{—p^) -C 1. Expanding the denominator in eq. ([T3|) around ^ = one obtains 
a sum over moments of the pion wave function: 

F*{p\ (p + qf) = ^T.C f'du u-^M ■ (15) 

From the definition (0) it is easy to see that these moments are given by vacuum-to-pion 
transition matrix elements involving increasing powers of the covariant derivative. For = 
[p + g)^, i.e. g = 0, only the lowest moment n = contributes in eq. ([T5|), and the form 
factor reduces to Fq/p^ which is the classical result. In contrast, if the photon virtualities 
differ strongly from each other, then many moments contribute to eq. ([T5|). In this case, the 
calculation of the form factor requires the knowledge of the shape of the pion wave function. 

Returning to the correlation function (^ one realizes that the same technique may be used 
to obtain a representation analogous to eq. ([13|) . The only new element in the correlation 
function (^ is the virtual heavy quark propagating between the points x and instead of the 
light quarks present in eq. (|12D. This gives rise to important differences which however do not 
change the formalism substantially. For the present discussion it is sufficient to stick to the 
approximation (§) and confine ourselves to the first term proportional to ttIc- The complete 
analysis of this expression and the calculation of further corrections will be carried out in the 
next section. Furthermore, writing in terms of invariant amplitudes: 

F^(p, g) = F{p\ (p + qf)q, + F^p", (p + qf)p^ , (16) 

we focus on the function F. Using the definition eq. ([T^ ) of the leading twist wave function 
and integrating over x and k one finds 

F(/, {p + qf) = mj. f /"^f ^ . (17) 

JO — [p + uqY 



Thus, the infinite series of matrix elements of local operators encountered before in eq. ([T0| ) 
is effectively replaced by an unknown wave function. The expression (0) is rather similar to 
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the one quoted in eq. ([T3|) for the TT^'~f*'y* form factor. Most noteworthy is the fact that the 
large- distance dynamics is described by one and the same pion wave function. This universal 
property is essential for the whole approach. 

Next we indicate how the relation ( p!7D can be turned into a sum rule for the coupling 
constant gD*Dn- The key idea is to write a hadronic representation of F by means of a double 
dispersion integral: 

m'jjmD*fDfD*gD*Dn 



F{p\{p + qy 



+ 



p^{si, S2)dSidS2 
P2{s2)dS2 

S2 - {p + qY ' 



+ 

p\{si)dsi 

Si 



The first term arises from the ground state contribution already indicated in eq. while 
the spectral function p'*(si,S2) is supposed to take into account higher resonances and con- 
tinuum states in the D* and D channels. The additional single dispersion integrals originate 
in subtractions which are generally necessary to make the double dispersion integral finite. 
Then, considering and {p + qY independent variables one can perform the usual Borel 
improvement in both channels. Applying the Borel operator 



lim. 



d 



dQ' 



(19) 



to eq. ([T8|) with respect to p and {p + q) \ we obtain 



F{MlMl) ^ BMiBMlF{p\{p + qf 



ml)mD*fDfD*gD*D7, 



TTlr 



+ 



p'^iSi, S2)dSidS2 , 



(20) 



where Ml and M| are the Borel parameters associated with p^ and {p + g)^, respectively. 
Note that contributions from heavier states are now exponentially suppressed by factors 

exp{ — ^'^ A,r^°*'° } desired, while the subtraction terms depending only on one of the vari- 



ables, or {p + qY, vanish. 

The same transformation has to be applied to the expression 
{p + uqY = (1 — u)p^ + u{p + qY, and use 

(/-I)! 



To this end we rewrite 



BmI^mI 



rat 



2\2-l^-m; 



(1 — u)p'^ — u{p + qYY 
where the Borel parameters and M| have been replaced by 

Ml ,,2 Ml Ml 



■""'Sin -no) 



Uo 



Mf + Mi 



2 ' 



M' 



Ml + M| 



Finally, equating the quark-gluon and the hadronic representations of F{Ml,Ml 
carding for a moment contributions of higher states, we end up with the sum rule 



rur. 



■ go'Dn = mj^ v9^(Mo)M^exp 



D* 



Ml 



+ 



m 



D 



Ml 



+ . 



(21) 

(22) 
and dis- 

(23) 
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The ellipses refer to higher-twist contributions which we discuss in detail later. Since and 
M| are expected to be quite similar in magnitude, the coupling constant gD*Dn is determined 
by the value of the pion wave function at m ~ 1/2, that is by the probability for the quark and 
the antiquark to carry equal momentum fractions in the pion This interesting feature is 
shared by the sum rules for many other important hadronic couplings involving the pion. 

The quantity (y9^(l/2) is considered to be a nonperturbative parameter, similar to quark 
and gluon condensates in the standard approach. It may be determined from suitable sum 
rules in which the phenomenological part is known experimentally. We use the value 



(^^(1/2) = 1.2 ±0.2 



(24) 



obtained in ref. |T^ . 

The dependence on the pion wave function disappears in the kinematical limit g — as can 
be seen from eq. (P^). This is just the limit where the correlation function (|^) can be treated 
in SDE. The condition g ~ is implicitly assumed in refs. |1I0| , [1^ where the correlation 
function is calculated using the external field method. This technique is equivalent to 
the soft-pion approximation used in refs. [pT] , p!3| as will become clear later. For comparison, 
we present the sum rule following from eqs. (p!7|) and (18) by putting g = 0, or equivalently 
(p + g)2 = Since is the only variable left, one now can only perform a single Borel 
transformation and, hence, the subtraction terms in the double dispersion relation (^) are no 
more eliminated. Moreover, the contributions to eq. (|T8|) of transitions from excited states to 



ground states are not suppressed after Borel transformation |g6|. This point will be explained 
in more detail later in Sect. 5 . In the approximation considered in eq. (EBf) one obtains 



■ go'Dn + M^A = rricUM^ exp 



m 



D* 



2M2 



+ 



m 



D 



2M2 



+ 



(25) 



where A is an unknown constant corresponding to the contributions of unwanted transitions 
and subtraction terms. 

From eqs. (p3[) and (25) one can clearly see the advantages and disadvantages of the 
two approaches. In the light-cone sum rule (^) the hadronic input is simple, whereas the 
theoretical expression involves a new universal nonperturbative parameter, namely (Pt,^{1/2). 
Just the opposite is the case for the sum rule (^) at g = 0. Here the QCD part is straightfor- 
ward, while the hadronic representation now involves an additional unknown quantity, which 
is non-universal and specific for this particular sum rule. A comparison of the results obtained 
in these two approaches should allow one to check the reliability and improve the accuracy of 
the predictions. 



3 Light-cone sum rule for Qd^Dti and Qb^b-k 

In this section we systematically derive the light-cone sum rule for the D*D7i and B*Btt 
couplings taking into account the two- and three-particle pion wave functions up to twist 4. 
First, we complete the calculation of the diagram Fig. la which represents the contribution 
from quark- antiquark wave functions. To this end we return to the expression (^. In the first 
matrix element we include the twist 4 corrections in addition to the leading twist term already 
given in eq. ([l^): 

(7r(g)|J(x)7;.75w(0)|0) = -iq,U C due'^'^^ (^pM + 9i{u) + 0{x^)) 

•J 



6 




(26) 



On the r.h.s. of this relation one sees the first few terms of the hght-cone expansion in x'^ 
of the matrix element on the l.h.s.. While ip-rt parametrizes the leading twist 2 contribution, 
gi and g2 are associated with twist 4 operators. In the second matrix element of eq. (||) we 
substitute 

Ifilu = -icTfiu + gfiu (27) 

and express the result in terms of the twist 3 wave functions and y^o- defined by the matrix 
elements 

(7^(g)|(i»^75u(0)|0)= ^\ " / due'^^^^ifpiu) 



and 



(7r(g)|d(a;)(T^i.75u(0)|0) = i{q^x^ - qyX 



rriu + rrid Jo 



6(m„ + md) Jo 

It should be noted that in eqs. ( ^^^91) the path-ordered gauge factors 



Pexp{igs / dax^A^{ax)} 





(28) 
(29) 

(30) 



appearing in between the quark fields and assuring gauge invariance, are not shown for brevity 
since they formally disappear in the light-cone gauge x^A^ = assumed throughout this paper. 



More details on these wave functions can be found in refs. 0, |2^ and in Appendix A. 

Collecting all terms , we obtain the following result for the invariant function F as defined 
in eq. (p!6|): 

F^^\p',{p + qf) = 



du 



ml — {p + uqY 
1 / p^ + m1 

6 I ml — {p + uqY 



mcULp^{u) + 



where 



G2(u) 



rriu + ma 

^ml{g^{u) + G2{u)) 

m^ — {p + uqY {m^ — {p + uqYY 



g2{v)dv 



2ug2{u) 



(31) 



The suffix (a) refers to the diagram in Fig. la which represents the leading twist term in the 
light-cone expansion of the c-quark propagator given in eq. (^. 

In addition, to the accuracy of eq. (|3TD we must also take into account higher twist terms 
in the propagator up to twist 4 which are numerous, in general. The complete expansion 
is given in ref. [27|. One has contributions from qOq, qGGq and qqqq nonlocal operators. 



G denoting the gluon field strength. Here we only consider operators with one gluon field, 
corresponding to quark- ant iquark-gluon components in the pion, and neglect components with 
two extra gluons, or with an additional qq pair. This is consistent with the approximation of 
the twist 4 two-particle wave functions derived in ref. and used here. Taking into account 
higher Fock-space components would demand corresponding modifications in the two-particle 
functions via the equations of motion. Formally, the neglect of the qGGq and qqqq terms can 
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be justified on the basis of an expansion in conformal spin |28|. In tliis approximation tlie 
c— quark propagator reads 



(0|T{c(x)c(0)}|0) = tS'',ix)-igs 



ikx I 



(27r 







1 ^ + rUc 
2(m2- A;2)2 



+ 



ml — A;2 



(32) 



wliere G^^ = G'^^^^ witli tr(A"A ) = 26 , and gs is tlie strong coupling constant. 

Substituting eq. ( ^2l) into eq. (|]) and using eq. (|l^) one obtains tlie contribution to the 
invariant function F represented by the diagram in Fig. lb: 



d Jv d X di) . I "7 , , 
t / , ■ , „ —{Tx\d[x)-i^ 



+ 



(27r)4(m2 - A;2) 



vXpG'^^{vx)'y\ 



ml — k'^2 



75w(0)|0) . 



With eq. (^Tf) and the identities 
and 

ifiiu(yp\ = {(^p\9up - (y^lpgu\) + i{9ti\9up - gppgux] 

one is led to the three-particle pion wave functions defined by 

{'K\d{x)gsGf,y{vx)aap'-f^u{Q)\Q) = 

= ihTT[{q,j.qagvp - quqagixp) - {.qi^qpgua - quqpgtxa)] j 'Dai(p3^{ai)e 



(33) 



(34) 



(35) 



(36) 



{Tc\d{x)-ff,-fr,gsGap{vx)u{0)\0) 

= u 



/ Xaqp,\ / xpqp 

qp ]^^o.\ 



qx 



qx 



Vaiip±{ai)e 



iqx{ai+va3) 



+ U^iqaXfs - qpx^) I Va, ^||(a,)e*^^(°^+''-3) 
qx ' 



(37) 



{■K\d{x)'jpgsGaf3{vx)u{0)\0) 



qp 9a,^ ] dPt^ 



qx 



qx 



Dai <f>±{ai)e' 



(38) 



where Gap = ^^aparG'^'^ and Vai = daida2da3S{l — ai — «2 — c^s)- While f3n{(yi) is a twist 3 
wave function, the remaining functions ip±, (p\\, ip± and ip\\ are all of twist 4. Substitution of 
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these expressions in eq. ([33|), finally, yields 



[ml - (p + («! + va2)qYY J 

In addition to Fig. lb there are further gluonic diagrams such as the ones depicted in 
Figs. Ic and Id. Note, however, that it is not necessary to take the diagram in Fig. Ic 
into account separately, since its contribution (to twist 4 accuracy) is already included in the 
two particle wave functions. In contrast, the two-loop perturbative corrections exemplified in 
Fig. Id should be included in a systematic way, but their calculation lies beyond the scope of 
this paper. 

Putting together eqs. ( pT]) and (p9| ) and applying the double Borel transformation pT| ) 



with respect to and {p + qY-, we end up with the following expression for the invariant 
amplitude F : 

F{Ml Ml) ^ F{M\ uo) = BmiBm?F{p\ {p + qf) = 

-""c of r /tt^^ / / \ 1 / \ 1 d(f„ ml \ 

~j4^^o32{uo) - -^^{giiuo) + G2{uo)) + 2/3^/3 [uo) + mj\, ^ . (40) 

Here, l9 and If involve the three-particle wave functions of twist 3 and 4, respectively : 



^Psniai, 1 - Uo, uq - ai) /"i""! v?3,r(ai, 1 - Qi - 0^3, Q3) 



If (no) = Hda, y-3.(ai, 1-^0, ^0-^1 j _ f 

Jo Uq — ai Ju 



Uq — Cii JuQ—ai 



(41) 



Jf (mo) = / dai —[2ifi^{ai) - ¥^||(a,) + 2^^{a,) - (^||(ai)] . (42) 

JO Juo-ai as 

The Borel parameters and uq are given in eq. (0). The above is the desired quark-gluon 
representation of the invariant amplitude F in the correlation function (^. 

The remaining task now is to match eq. ( ^OD with the corresponding hadronic representa- 
tion (1^) and to extract the coupling Qwd-k- As usual, invoking duahty, we assume that above 
certain thresholds in si and S2 the double spectral density p'*(si,S2) associated with higher 
resonances and continuum states coincides with the spectral density derived from the diagrams 
in Fig. 1. The procedure is explained in detail in Appendix B. For Ml = M| = 2M^ and 
Mo = |, and for standard polynomial wave functions, the effect of the continuum subtraction is 
remarkably simple |T^, . It amounts to the following replacement of the exponential factor 



multiplying the twist 2 and 3 terms proportional to in eq. (^Ol): 

(43) 



„2 



So being the threshold parameter defined in eq. ( |B10| ). The higher twist terms which are 
suppressed in eq. (EO) by inverse powers of with respect to leading one remain unaffected. 
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With eqs. {^^, (|40|) and (|43|) it is then easy to derive the following QCD sum rule for the 
DD*TX coupling: 



fofD-gD-D-K 



-e 2M 



2 



M [e m2 — e 



rric \ 3 6 du I rricU 



1 (iv^^ 



2/3. 



+ e 



V9^(uo) + 2uog2{uo) 



4ml 
M2 



((7iK)+G'2K)) + /fK 



where 



m: 



-2(gg) 



"0=1/2 



(44) 



(45) 



friu + nid /2 

In eq. (^) we have used the familiar PCAC relation between m^, , the quark masses and 
the quark condensate density {qq). Note that the twist 2 and 3 and the twist 4 wave functions 
have different dimensions (see Appendix A). G-parity implies (72(1/2) = dip^/ du{l/2) = 0, so 
that these terms vanish in the sum rule (|4^. 

For completeness, we also repeat the standard two-point sum rules for the decay constants 
fo and fo' ■ 



flrn 



D 



mi 



so '"j- 

dse M'-i 



m 



2\2 



1 + 



^0 
2M2 



and 



j:2 2 



57r^ Jm] 



'0 

dse 



mi , 
2M2- 



2 2 ^ 
m^m^ 

4M4 



(46) 



(47) 



where mg = {qCapG^^ q) / {qq) is a conventional parametrization for the quark-gluon conden- 
sate. In the above, we have omitted numerically insignificant contributions of the gluon and 
four-quark condensates. For consistency, we do not take into account perturbative 0{as) 
corrections to these sum rules, since they are also not included in the sum rule 



4 Numerical analysis 



The principal nonperturbative input in the sum rule (H3) are pion wave functions on the light- 
cone. In ref. a theoretical framework has been developed to study these functions. In 
particular, it has been shown that the wave functions can be expanded in terms of matrix 
elements of conformal operators which in leading logarithmic approximation do not mix under 
renormalization. For example, for the leading twist pion wave function one finds an expansion 
in Gegenbauer polynomials. 



99^(m, /i) = 6m(1 -u) 1 + a2(/i)C2 



2u-l) + aAii)cT{2u - 1) + 



(48) 



where all the nonperturbative information is included in the set of multiplicatively renormal- 
izable coefficients a„, n = 2,4, . . .. The corresponding anomalous dimensions are such that 
the coefficients a„ vanish for fi —>■ 00, and the wave function is uniquely determined by the 
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first term in the expansion. Therefore, this term is called asymptotic wave function. Similar 
expansions also exist for the wave functions of nonleading twist . 

For practical applications it is important that the expansion in conformal spin converges 
sufficiently fast. How fast the wave functions approach their well-known asymptotic form is 
still under debate. However, there are indications that the nonasymptotic deviations have 



been overestimated previously. Corrections to the asymptotic expressions in next-to-leading 
(and in some cases also next-to-next-to-leading order) in conformal spin are known for all 
wave functions which appear in the sum rule (@). For details, we must refer the reader to the 
original literature [1^, |28[ . 

In our numerical analysis we use the set of nonleading twist wave functions proposed in 
ref. The explicit expressions and the values of the various parameters are collected 

in Appendix A. Furthermore, we take = 132 MeV, fiT,{lGeV) = 1.65, corresponding to 
(qq) = —(243 MeV)^, mg = 0.8 GeV^ and, in the charmed meson channels, rric = 1.3 GeV, 
So = 6 GeV^, rriD = 1-87 GeV and niD* = 2.01 GeV. The same parameters are also used to 
determine the decay constants fo and fo* from the sum rules (^) and (^Tj). One obtains 

/d = 170 ± lOMeV, fo" = 240 ± 20MeV . (49) 

The uncertainty quoted characterizes the variation with the Borel parameter in the interval 
1 GeV2 < M2 < 2 GeV2. 

Having fixed the input parameters, one must find the range of values of for which 
the sum rule (|4^ ) is reliable. The lowest possible value of is usually determined by the 
requirement that the terms proportional to the highest inverse power of the Borel parameter 
stay reasonably small. The upper limit is determined by demanding that the continuum 
contribution does not get too large. In the D*D7t sum rule we take the interval 2 GeV^ < 
< 4 GeV^. In this interval the twist 4 term proportional to does not exceed 5%. 
Simultaneously, higher states contribute less than 30%. The dependence of the r.h.s. of eq. 
(^^ on the Borel parameter is shown in Fig. 2a. As can be seen, in the fiducial range of 
given above, the sum rule is quite stable. From Fig. 2a one can directly read off the prediction 

fofD-QD^D. = 0.51 ± 0.05 GeV' . (50) 

Dividing this product of couplings by the decay constants (^) finally yields for the D*D7r 
coupling constant: 

gn*Dn = 12.5 ± 1.0 , (51) 
where the error is understood to indicate the range of values corresponding to the correlated 



variation of the results (|49D and (|50D within the fiducial intervals. While the combination of 
couplings ( |50D is not affected by uncertainties in the decay constants, it is more sensitive to 
the charm quark mass and threshold sq than the coupling gD*D-K itself. 

A few comments are in order. The twist 3 terms contribute to the sum rule at the level 
of (50 -^ 60) % and are therefore as important as the twist 2 contribution. On the other 
hand, the impact of twist 4 is small amounting to about 5%. Two sources of uncertainties not 



included in eq. ( ^0]) are the nonasymptotic corrections to the leading twist wave function (^^^ 
and to the three-particle (twist 3) wave function (p^T,. The latter in turn induce modifications 
in the two-particle (twist 3) wave functions ipp and ip^, apart from the corrections generated 
by the asymptotic (p^^^ wave function. In order to estimate the sensitivity of our results to the 
nonasymptotic effects in yj^ and ^93^ we drop these corrections altogether and recalculate the 
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product fDfD'QD'D-K- Remarkably, the result changes by only 5%. One can thus be confident 
that the total uncertainty in eq. (|50D does not exceed 20%. 

We should emphasize that the above prediction can be directly tested experimentally in 
the decay D* — > Dti. With the value of gD*D-K given in eq. (0) one predicts the decay width 

rm*+->DV) = ^^^IgP = 32±5keV. (52) 

247rm|)* ^ ^ 

Predictions for other charge combinations are easily obtained from eq. ( |5^ ) taking into account 
the isospin relations @) as well as small differences in the phase space volumes: 

r(L)*+ ^ D°7r+) = 2- l.ir(D*+ ^D+7r°) = 2-0.72r(L)*° ^D°7r°) . (53) 

The current experimental upper limit 

T{D*+ D°7r+) < 89 keV (54) 



is obtained by combining the limit Ttot{D*^) < 131 keV [^] with the branching ratio 
BR{D*+ D^n+) = (68.1 ± 1.0 ± 1.3)% [0. Our prediction is well below this upper limit. 

The sum rule for qd* d-,^ given in eq. ( PD is easily converted into a sum rule for the coupling 
93* Bn = 9b*ob~-k+ by replacing c with 6, D with 5, and D* with B* . The corresponding 
parameters are = 5.279 GeV, niB* = 5.325 GeV, rrih = 4.7 GeV, and sq = 35 GeV^. In 
addition, one has to evolute the wave function parameters to a higher scale /if, (see Appendix 



A). With these changes the two-point sum rules ( ^61) and (^Tf ) yield 



fB = 140 MeV , fs' = 160 MeV (55) 

with negligible uncertainties due to variation of M^. Using the same criteria as for the go'Dn 
sum rule the fiducial range in turns out to be 6 GeV^ < < 12 GeV^. The stability of 
the sum rule for the product of couplings fBfB*gB*Bn is illustrated in Fig. 2b. We obtain 

fBfB'9B'B. = 0.64 ± 0.06 GeV' , (56) 

and with eq. (|55|) 

gB'B. = 29 ± 3 . (57) 

The hierarchy of various twists as well as the uncertainty due to the nonasymptotic corrections 
are found to be similar as in the case of gD*D-K ■ Contrary to the latter, the coupling constant 
Qb'Btt cannot be measured directly, since the corresponding decay B* Bti is kinematically 
forbidden. However, the B*Bti on-shell vertex is of a great importance for understanding of 
the behaviour of heavy-to-light form factors as will be discussed in Sect. 6. 



5 Sum rules from the short-distance expansion 

With the results of Sect. 4 at hand we are now also in a position to study in more detail the 



soft-pion limit ( pS]) of the sum rule (^) which is obtained from the correlation function (|^) at 
= {p + qY or, equivalently, for g 0. As discussed in Sect. 2, in this limit one can apply 
a short-distance expansion in terms of local operators with increasing dimensions in contrast 
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to the light-cone expansion involving nonlocal operators with increasing twist. In technical 
terms, at g — only the lowest moments of the wave functions contribute. Thus the integrals 
reduce to overall normalization factors. The explicit expression for the invariant function F 
in this limit is directly obtained from eqs. ( |3TD and (p9|): 



ml — 



1 + 



2 + 



mt. 



ml — 



+ 



105^ 



9(m^ — p"^) 



mt. 



ml — p^ 



(58) 



where the parameter 5 is specified in Appendix A. After Borel transformation in p one has 



1 f u-wmr 10 5m?(5^ 



1 + + 

3m, M2 



£ + _52 

9 7 9M4 



(59) 



As indicated in eq. (pq), the price for simplifying the QCD representation of the correlation 
function is a more complicated hadronic representation. This in turn makes it more difficult 
to extract the ground state contribution containing the D*Dn coupling. For illustration, we 
consider the contribution in the dispersion representation of the correlation function (^ from 
the transition of a given excited state in the D*-channel with mass > m^* to the ground 
state Z)-meson at g ^ 0. This contribution is proportional to 



(ml — p'^){m?j^ — p2) 



and, after Borel transformation, to 



(60) 



ml — m\) 



/ ™2 2 



(61) 



Similar expressions hold for the ground state transition D* 
mo = rriD*, one has a double pole instead of eq. (|5I1|) and 



D with m^, = m£,*. In the limit 



1 

M2' 



(62) 



instead of eq. (^). Clearly, the contribution (|6lD is not exponentially suppressed relative to 
the ground state contribution (^), and can therefore not be subtracted assuming duality. On 
the other hand, transitions involving excited states in both the initial and the final states are 
suppressed by Borel transformation with respect to the ground state transitions and cause no 
problems. Schematically, the complete hadronic part of the invariant amplitude F{M'^) can 
be written as follows: 



1 f ml,mD*fDfD' 



M2 \' 



go'Dn + AM^ e 



4 



+ C-e 



(63) 



where the constant A incorporates all unsuppressed contributions of the type (|6ll) , while the 
term proportional to C contains all exponentially suppressed contributions. 

To get rid of the contaminating term AM"^ we follow ref. and apply the operator 



1-M^ 



d ' 



2 I 2 
2M^ 



(64) 
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to both representations ([59|) and (|6^) . Now it is possible to subtract the continuum contribu- 



tion by the same substitution (^31) which we have already employed to get the light-cone sum 
rule (|4^). In this way we obtain a new sum rule for gD*D-K'- 



M'il - 1 + ^ + + ^8 



3m, y 3 9 9M2 



(65) 



For the same input parameters and the range of leading to the prediction (^Tj) this sum 
rule yields 

= 11±2. (66) 

Similarly, replacing in eq. (|65D the charmed meson parameters by the corresponding quantities 
in the beauty channel, one finds 

gB*B. = 28 ± 6 . (67) 

As compared with the predictions (^) and ( |^ the uncertainties are larger by a factor of 
two due to the worse stability of the sum rule ( |65D against variation of M^. The agreement 
of the results indicates selfconsistency of the sum rule approach and gives support to the 
approximations used in the pion wave functions. 

Furthermore, one can show that the sum rule (|65D is actually equivalent to the sum rule 



obtained in ref. |10] using external field techniques. Indeed, applying the usual reduction 



formalism to the pion, one can rewrite the correlation function (^) in the following form: 

FAV. q) = iq' - O J d'x d% e'(P^-^y\0 \ T{J(x)7,c(x), My), 0(0)75^(0)} I 0) , (68) 
where 07r(y) is the interpolating pion field. According to PCAC 

jfiiy) = ^{y)l iJ.l^div) being the axial current. Subsituting eq. (^) in eq. ( |68D and integrating 
by parts, one gets □ 

2 2 

— 777 ? 

F,{v,q) = ^^^^T,,{p,q)q'' = --T,,{p,q)q^ , (70) 

with 

T,,{p,q) = /rf'x Ae^^"\0|T{J(a;)7^c(x),j,^(y),c(0)75n(0)}|0). (71) 

Instead of dealing with the three-point correlation function T^p{p,q) directly, it is more 
convenient to consider the following two-point correlation function in the constant external 
axial field A^: 

T^{V, q)= I d'x e'P^iO \ T{d{xh^c{x), c(0)75«(0)} | 0)^ . (72) 



^In addition, one obtains two-point correlation functions because of contact terms. These do not lead to 
double poles in in the relevant dispersion relation and are therefore eliminated by applying the differentiation 



operator (|64p. 
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It is assumed that a term A^f^j^ corresponding to the interaction of the external field with 
the light quarks is added to the QCD Lagrangian. To first order in the external field, this 
correlation function is given by 

T>,g)=T^,(p,g)A^^ 

with Tf^p{p,q) as defined in eq. ([71|). In the above sense, the two-point correlation function 
(|72|) in the constant axial field, is equivalent to the three-point function (| 7ll) and, via the 
PCAC relation (|69|), also to the two-point correlation function (^) at g — a. Consequently, 
the sum rules obtained in refs. and |jl3| should coincide with each other and with the sum 



rule ( [65| ) derived in this paper. 

In particular, the expression ( ^9]) for F(M^) can be compared with the result given in eq. 
(19) of ref. [jlO| and in eq. (2.15) of ref. |13[, after normalization and kinematical structures 
are adjusted properly. In refs. [|10|, the correlation function (^ is separated as follows: 

= Aq^ + B{2p^ + q^) , (73) 

and the sum rule is obtained by evaluating of the invariant function A for g ^ 0. In terms of 
the invariant amplitudes defined in eq. ([16D one has 



A = F-|^, 5 = ^. (74) 

Hence, for comparison we need also the second invariant function F, which we have not 
discussed so far, but which can be calculated along the same fines. Adding this contribution 
to eq. (|59D , we have checked that our result for the invariant amplitude A coincides with 
the one presented in ref. |T3|, apart from terms proportional to ttt-q which, being associated 
with twist 5 contributions in the light-cone sum rules, are neglected in our approximation. 
Numerically, this terms are not important. On the other hand, we disagree with ref. [|T^] in 
the non- leading terms proportional to 5^. 

Although it is legitimate to use different Lorentz decompositions of the correlation function 
(0) in order to derive the desired sum rule, we think that the choice adopted in the present 
paper is more adequate for the following reason. Since the vector current g7^c is not conserved, 
it not only couples to = 1~ vector mesons, but also to = 0"*" scalar mesons (-Do)- The 
corresponding transition matrix element is proportional to the momentum p^: 

(0 I g7^c \ Do) = fDomDoPfi ■ (75) 

The mass of the ground state Dq meson is expected to be in the vicinity of 2400 MeV which is 
not far from the mass of the D* and below the accepted continuum threshold in the D* channel. 
For this reason, the Dq contribution should be added to the sum rule. Unfortunately, this 
introduces additional uncertainties in the hadronic representation as is the case, for example. 



in ref. |jT3[. In contrast, our sum rules based on the invariant function F defined in eq. ([T6| ) 



are not affected by scalar meson contributions, which is a clear advantage. 



A calculation rather similar to ref. 13 1, but with particular emphasis on the heavy quark 



limit, has been carried out earlier in ref. |11|. Very recently, another calculation in the heavy 



^ As a side remark, the light-cone approach leading to eq. ( [44| ) corresponds to a calculation in the background 
of a variable external axial field |T^, . 
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quark limit appeared in ref. ||T2| using the external field technique. Unfortunately, in this 
paper a wrong expression for the induced quark condensate in the external field is used, as 
can be seen by consulting refs. The error can be traced back to the equation of 

motion for the quark field in presence of the external field which is modified from il^q = to 
ilj^q =4^^q. By this modification the axial current insertions into the vacuum quark legs are 
properly taken into account. The numerical comparison of these different calculations is left 
for the concluding section. 



6 Pole model for D ^ tt and B ^ tt form factors and 
QCD sum rules 

The couplings Qo'Dtt and gB*B-K fix the normalization of the form factors of the heavy-to-light 
transitions D ^ it and B tt, respectively, in the pole model description [5,6]. This model 
is based on the vector dominance idea suggesting a momentum dependence dominated by the 
D* and B* poles, respectively. More definitely, the form factor f^ip^) defined by the matrix 
element 

(7r(g) I d^^c I D{p + q)) = 2fUp')q, + UUp") " /d(p')K (76) 
is predicted to be given by 

An analogous expression holds for the form factor fsip^)- 

It is difficult to justify the pole model from first principles. Generally, it is believed that 
the vector dominance approximation is valid at zero recoil, that is at mj^. Arguments 

based on heavy quark symmetry suggest a somewhat larger region of validity characterized by 
(m|,— p^)/mc ~ O(lGeV). However, there are no convincing arguments in favour of this model 
to be valid also at small values of which are most interesting from a practical point of view. 
Therefore, the finding 0, Q that the pole behaviour is consistent with the dependence at 
p'^ predicted by sum rules, is very remarkable. Meanwhile, this claim has been confirmed 
by independent calculations within the framework of the light-cone sum rules Q. 

In this section we want to demonstrate that not only the shape but also the absolute 
normalization of the above form factors appears to be comparable with the pole model de- 
scription. This assertion is non-trivial, since contributions of several low-lying resonances in 
the D* or B* channel could still mimic the p"^ dependence of a single pole, but the relation 
to the coupling gD*D-K or gB*B-K should then be lost 0. However, despite of the overall agree- 
ment in the mass range of D and B mesons, there is a clear disagreement on the asymptotic 
dependence of the form factors on the heavy mass. The QCD sum rules on the light-cone 
provide a unique framework to examine these issues, since both the form factors fDsiv^) 
''^1b~ ^ 0(1 GeV^) and the couplings Qo'D-k and gB*B-K can be calculated from the same 
correlation function (|]) using the same technique. In addition, contrary to conventional sum 
rules [0, this approach leads to consistent results in the heavy quark limit ||21|| . 



The detailed derivation of the light-cone sum rules for the D — > vr and B ^ n form factors 
is discussed in ref. (see also refs. [19-21]). Here we just mention that the sum rule for 
/d(p^) is obtained by matching the expressions (0) and (|39D for the invariant amphtude 
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{p + g)^) in terms of the pion wave functions with the hadronic representation 



2mlfnfMp') 

^c(^D ~ {P + 



+ 



°° p^{p'^,s)ds 
s - {p + g)2 



(7J 



In the above, the pole term is due to the ground state in the heavy channel, while the excited 
and continuum states are taken into account by the dispersion integral above the threshold 
So- Invoking duality, the latter contributions are cancelled against the corresponding pieces in 
eqs. ([31| ) and ([39|). After Borel transformation in the variable {p + g)^, the resulting sum rule 
takes the form 



"1 du 
— exp 



m 



D 



— — u) 



M2 



udu 



VaiQ{ai + uas — A) 



exp 



m 



D 



M2 



$2(M,M2,p2) 

- p^il — ai — ua^) 
{ai +ua^)M'^ 



(79) 



<l>3(M,M2,p2 



where 



$2 



V?vr(M) + 



ml + p^ ' 
mM2 



Amlg^{u) 2G2{u) 



uA/P 



1 + 



+ p 



^0) 



2/3. 



2 2 



(ai + Ma3)M2 
2v2±(ai) - V'ii(ai) + 2(^_L(ai) - <^||(ai) 



and A = (m^ — p'^)/ [sq — p^). The notation is as in eq. (|4^). Improving the approximation 
given in ref. we have added the contributions of three-particle wave functions of twist 4. 
The analogous sum rule for the B ^ n form factor follows from the above by replacing c ^ b 
and D ^ B, and by rescaling /i^ and the wave function parameters from /i^ to /if,. 

The maximum momentum transfer p"^ at which these sum rules are applicable is estimated 
to be about 15 GeV^ for B mesons, and 1 GeV^ for D mesons. For numerical evaluation we 
use the approximations of the wave functions given in Appendix A. We emphasize that the 
input here is exactly the same as in the calculation of the couplings gD*D-K and gB*Bn ■ The 
form factor f^ip^) resulting from the sum rule (pOf ) is plotted in Fig. 3a, together with the 
corresponding prediction (|77D of the pole model. We see that in the region of overlap both 
calculations approximately agree with each other. To a lesser extent, this also applies to the 
form factor fsip^) illustrated in Fig. 3b @. Quantitatively, at = we find 



/^(0)5ii = 0.66, /+(0)PM = 0.75 , 



(82) 



and 



/l(0)5i?, = 0.29, /+(0)pAf = 0.44 . (83) 

Thus, in the regions tUq — p^ > 0(1 GeV^) with Q = c and b, respectively, the numerical 
agreement between the light-cone sum rule and the pole model is better than 15% for f^, but 
only within 50% for f^. 



*The dependence of eq. (|8^) on the Borel parameter is weak For definiteness, we take here — 4 
GeV^ for the D ^ tt form factor and A'P = 10 GeV^ for the B tt form factor. 
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At this point, we must add a word of caution concerning the apphcabihty of the pole model 
too far away from the zero recoil point, in particular at = 0. The two descriptions differ 
markedly in the asymptotic dependence of the form factors on the heavy mass. Focusing on 
B mesons and using the familiar scaling laws 

fBy/rn^ = /b, fs^V^ = fs- , (84) 

and 

9b*b-k = —r- ■ 9 , (85) 

Jtt 

which are expected to be valid at rrih — > oo modulo logarithmic corrections, one obtains 

fUO)pM ~ l/V^ , (86) 
whereas the light-cone sum rule (79) yields [|l^] 

fUO)sR ~ l/mf . (87) 

This result rests on the QCD prediction of the behaviour of the leading twist pion wave 
function near the end point, that is (pniu) ~1— MatM^l. It should be noted that the 
contribution estimated by the sum rules corresponds to the so-called Feynman mechanism. In 
the case of heavy-to-light transitions it leads to the same asymptotic behaviour as the hard 
rescattering mechanism |T^, Recently it has been shown that the power behaviour 



(|87D of hard rescattering is not modified by the Sudakov type double logarithmic corrections. 
We believe that the light-cone sum rules reproduce the true asymptotic behaviour, although a 
rigorous proof in QCD is still lacking. On the other hand, we see no theoretical justification for 
extrapolating the pole model to the region p"^ = 0. The solution suggested by Fig. 3 is to match 
the two descriptions in the region of intermediate momentum transfer ^ rriq — O(lGeV^). 



Referring for a detailed discussion to refs. [21] and |3H| we want to emphasize that the light 



cone sum rules seem to be generally consistent with the heavy quark expansion. In particular, 
the light-cone sum rule (|4^) correctly reproduces the heavy quark mass dependence of the 
coupling gB*Bw Fitting our predictions for Qb'Btt and gD*Dn to the form 



2m B ^ 
Qb-Btt = — • g 

Jtt 



A 

rriB 



and the analogous expression for gD*DiT, we find for the coupling g and the strength A of the 
l/rnq correction: 

^ = 0.32 ±0.02, A = (0.7 ±0.1) Gel/ . (89) 

Furthermore, we are able to make a numerical prediction for the theoretically interesting 
ratio 

gB^B.fB*V^ ^ Q g2 . (90) 

go'DnlD'y/rnB 

This ratio is unity in the heavy quark limit and is shown to be subject to l/rriQ corrections 
only in the next-to- leading order . Our result (^Dj) is nicely consistent with this expectation. 
The deviation from unity also agrees in magnitude with the estimate in ref. [0, but has a 



different sign. This is due to a sizeable difference in the ratio Jb* / fo* ■ While the values of 
the decay constants given in eqs. (|i7| ) and (^) yield 

^^^^ = 1.12 (91) 
fD*^/rnB 



18 



in agreement with the expectation quoted in ref. [Q, the latter ratio turns out to be larger by 
30% if calculated from fs* and fa* as assumed in ref. |T3 . 



7 Summary and conclusions 

We have presented a comprehensive analysis of the pion couplings to heavy mesons in the 
framework of QCD sum rules. The main new result of this paper is the light-cone sum rule 
providing the numerical estimates for go'Dn and gs's-K given in eqs. ( [5T| ) and (|57D, 
respectively. The decay width V{D* Dn) predicted in eq. ( ^2|) turns out to be three times 
smaller than the present experimental upper limit. We have compared our results to earlier 
QCD sum rule calculations [10-13], and resolved the existing discrepancies. 

A rather complete compilation of estimates i on the pion couplings to heavy mesons is 
given in Table 1. In the first row we show predictions on the reduced coupling g defined in eq. 
(1851). As one can see, the values obtained by combining the nonrelativistic constituent quark 
model with PCAC [H, |3^, ^ are roughly two times larger than the values favoured by our 



sum rule. However, more recent analyses [^, combining chiral HQET with experimental 



constraints on D* decays tend to give somewhat smaller values of g. Moreover, another recent 



calculation based on the extended Nambu-Jona-Lasinio model and chiral HQET is in 
perfect agreement with our estimate. 

The next two rows list the estimates of the couplings gB*B-K and go'Dn- These predictions 
are even wider spread. Quark models [^, ^ seem to give the strongest couplings, whereas 
SU(4) symmetry |44| and the reggeon quark-gluon string model [4^ predict a relatively small 



coupling. Two comments are in order concerning the analysis of ref. Firstly, these 

predictions are consistently lower than ours. There are several reasons for that: the different 
Lorentz decompositions (|l^) and ([73|) of the correlation function (||), the differences between 
the sum rule (Q) and the soft-pion limit ( |65D of it, the different regions of the Borel parameter 
M^, and the different values used for the decay constants //)(.) and fsi*)- In fact, as can be 
seen in Fig. 2, the couplings shrink with M^. However, given the reliability criteria, generally 
accepted for sum rules, we see no possibility to shift to larger values beyond the regions 



considered in this paper, in contrast to ref. . Secondly, we find it inconsistent to include the 
perturbative gluon correction in the estimates of fD,D* and fB,B* , since they are not included 
in the sum rule for the combination of couplings fofD^gwD-K and fsfB'gB'B-K- At least, we 
see no convincing argument in favour of such a procedure. For these two reasons we believe 
that the couplings are underestimated in ref. [l^ . 



For convenience and direct comparison with future measurements the decay width 
V{D*^ D^TT^) as calculated from gD*D-K oi g is quoted in the last row of Table 1. The 
widths in the channels D*^ D^n^ and D*'^ D^n^ are related to the above by coefficients 
which can be read off from eq. (|53D . Note that in contrast to the evaluation of r{D* Dn) 
from gD*Dn in this paper and in ref. |TB| the estimates in refs. B^, B^] using the reduced 



coupling g do not include 1/m corrections. However, the latter are important as can be seen 
from eq. (^). 

In addition, we have examined the pole model for the B ^ n and D ^ n form factors. 
Using our results on the gB'B-w and go'D-w coupling constants, we have found approximate 
numerical agreement between the pole model description and the direct sum rule calculation. 



^We have not included the results of ref. (lO| since to our knowledge this analysis is being reconsidered 
|. The result of ref. u% is omitted for reasons explained in Sect. 5. 
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However, the dependence on the heavy quark mass is found to be completely different in the 
region of small momentum transfers. We have argued in favour of the sum rule approach. 
Moreover, writing a heavy quark expansion for the couplings gB*B-K and Qd^d-k we have deter- 
mined the expansion coefficients, in particular, the leading 1/m correction. 

Last but not least, we have discussed in some detail the theoretical foundations and ad- 
vantages of the light-cone sum rules, complementing the work of refs. [14-19]. We believe that 
this approach is especially suitable for the study of heavy-to-light decay form factors, and 
coupling constants of the type considered in this paper. Further obvious applications include 
the radiative decays D* —>■ Dj and B* B'y. Since the photon wave functions are expected 
to deviate less from their asymptotic forms than the pion wave functions [0, these decays 
should provide a rather conclusive consistency check of the light-cone approach. 
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Appendix A 

For convenience, we collect here the explicit expressions for the pion wave functions used in 
our numerical calculations and specify the values of the parameters involved. 



For the leading twist two wave function we take [|T8| 

if^iu, /i) = 6m(1 -u)[l + a2ifx)C2^'^i2u - 1) + a4(Ai)Cf ^(2m - 1)] (Al) 
with the Gegenbauer polynomials 

C',/\2u-l) = ^[5{2u-lf-l] , 

Cf\2u - 1) = ^[21 (2m - 1)^ - U{2u - 1)2 + 1] , (A2) 
8 

and the coefficients a2 = |, 04 = 0.43 corresponding to the normalization point /i = 0.5 GeV. 
In the present applications the appropriate scales are set by the typical virtuality of the heavy 
quark. We choose 



= ^mfj -m2 - 1.3GeV, /ife = - ~ 2.4 GeV . (A3) 

Renormalization group evolution of the coefficients 02 and 04 to these higher scales yields 

a2(/Uc) = 0.41, a^ifJ^c) = 0.23 , 

a^ifib) = 0.35, a^ifib) = 0.18 . (A4) 

We stress that the value of <y9^ at m = 1/2 varies by only 2% when the scale /z is increased 
from 0.5 GeV to 2.4 GeV. Obviously, one can neglect this effect given the 15% uncertainty in 
the value of (p^iu = 1/2, fi = 0.5 GeV) quoted in eq. (p^). 
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According to the analysis in refs. [jl8|, g§] the set of wave functions of twist three is 
uniquely specified by the choice of the three-particle wave function (p^^^. Taking into account 
the contributions to ip^.-^ up to next-to-next-to-leading order in conformal spin, one has 



360a 



10203(1 +u;i,o^ (Tag - 3) 



1-^2,0(2 — 4aia2 — 80:3 + 803) + u;i,i(3Q;ia2 — 2^3 + 303 



This implies for the two-particle wave functions of twist three p8 

ifpiu) ■ 



1 + B2-{3{u - uf - 1) + 54-(35(m - - 30(m - uf + 3) 
2 8 



and 



6uu 



1 + C2^(5(^i - uY - 1) + Ci^i21{u - - U{u - u)^ + 1) 
2 8 



(A5) 



(A6) 



(AT) 



where 



B2 

C2 
C4 



30R, 
3 

^(5 - ^^1,0) , 



with 



R 



(A8) 



(A9) 



The coefficients /g^^ and Ui^k have been determined at the normalization point = 1 GeV from 
QCD sum rules []TB|: 



fsn = 0.0035 GeV^ tui,o = -2.88 , cj2,o = 10.5 , cui,i = 



(AlO) 



-1.05 



After renormalization ||23] to the relevant scales ( |A3D , we get 

/3.(/ie) = 0.0032 GeV^, a;i,o(/ic) = -2.63, a;2,o(/ic) = 9.62, ^1,1 (/ij ^ 

fsAl^b) = 0.0026 GeV^, luiM = -2.18, cu2,o(/ife) = 8.12, uji,iifib) = -2.59 

The corresponding numerical values of the coefficients ( |A^ ) follows: 



(All) 



52(/Xe) = 0.41, S4(/i,) = 0.925, C2 (/ic) = 0.087, C4(/ic) = 0.054 , 
B2{nb) = 0.29, 54(^6) = 0.58, C2{iJb) = 0.059, C4(/ib) = 0.034. 



(A12) 



In addition, the running of light-quark masses induces a scale dependence of the parameter 



/i^(l GeV) = 1.65 GeV, /i^(/ic) = 1.76 GeV, fi^{fib) = 2.02 GeV . 



(A13) 
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The wave functions of twist four are more numerous. The complete set given in ref. ^8 



(see also ref. I^Sj) includes four three-particle wave functions. However, in leading and next- 
to-leading order in conformal spin, these are specified by only two parameters: 

(fwiai) = 1205^e(ai — a2)aia2a3 , 

(^x(a^) = 3052^2(1 - as) [i + 2£(1 - 203)] , 

^\\{ai) = -1206^a,a2as[^ + e{l-3a3)] . (A14) 

The two-particle twist 4 wave functions are related to these by equations of motion. To the 
above order in conformal spin they involve no new parameters and are given by 

5 1 6 

g^(u) = -5^u\^ + -e6^[uu{2 + 13uu) + 10u^\nu{2-3u + -u^) 

6 

+ 10m^ lnM(2-3M + -M^)], 

5 

G^iu) = ^dW^ . (A15) 
One of the parameters is defined by the matrix element 

{7r\gjG^^ru\0) = tS^Uq^ . (A16) 
The QCD sum rule estimate of ref. yields 6^ = 0.2 GeV^ at = 1 GeV. The remaining 



parameter is associated with the deviation of twist four wave functions from their asymptotic 



form. At /i = 1 GeV it takes the value e = 0.5 p8[. Renormalization to the relevant scales 
Q gives 

6\fic) = 0.19, =0.45, 

S^ifib) = 0.17, = 0.36 . (A17) 

This completes the description of the pion wave functions, as far as it is needed for the 
apphcations in this paper. 

Appendix B 

Here we derive the substitution (^) used in the sum rule (0) in order to subtract the con- 
tinuum contribution. For this purpose we have to write the invariant amplitude F given by 
eqs. (^) and in the form of a double dispersion integral: 



/m2 Si - p2 7^2 {p + qY 

Focusing first on the leading contribution ([T7|): 



F{p\{p + qf)=mjJ ^—pL^L^ = mjJ , . ^ ^ r, (B2) 

Jo m'^^ — [p + uqY Jo — + gj^M — p"^(l — -u) 
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and changing variable from u to (m^ — p'^)/{s — p"^), one obtains 

{p + qf) = m J. r . f^-(^y^ . (B3) 

[s - [p + qy){s - p^) 

In general, the wave function ipn{u) can be expressed as a power series in (1 — u): 

k k \ ^ p J 

Substituting this representation into eq. (p3| ) and introducing formally two variables Si and 
S2 instead of s, it is easy to rewrite this expression in the form ( |B1| ) with the double spectral 
density 

p^^^is,, s,) = mj^ E f^rnj^'^ " - . (B5) 

The validity of eq. ( P5|) can easily be checked by direct calculation. The above derivation 
may seem tricky. However, there is a convenient general method to find the appropriate 
double spectral densities. One takes the Borel transformed invariant amplitude F(M^,M|) 
and performs two more Borel transformations in the variables ti = and T2 = 1/M|, to 

get 

B^,B,,F{l/n, I/T2) = p«^^(lM, 1/^2) . (B6) 

Details and useful examples can be found in ref. [^] . 

To proceed, we apply a double Borel transformation to the dispersion integral with 
pQCD gjygj^ gq (^p5[): 

B,,.B,riF = mj^ E r r ^^^^Asi - m^f5^'\s, - s,)e~^^/''le~^^l'''^ . (B7) 

^ ^ V Jml Jml 1 (fc + 1) 

Introducing again new variables s = S1 + S2 and w = si/s we can use the 5-function to evaluate 
the integral over v. The result is 



k 



m 

V — 



2~ 



k 



( svMl + s{l-v)Ml 

X GXp I 



(B8) 

J v=l/2 



At Mf = M| = 2M^ the f-dependence of the exponent disappears and the differentiation of 
the bracket gives a factor k\. We then get 

F(M2) = mj^ Y^^iT ds = m,U^^{l/2)M\''^ , (B9) 

which is the leading contribution in the sum rule (0). For arbitrary values of Mf and M| a 
similar expression is obtained, with the argument of the wave function and the Borel parameter 
in eq. ( [B9D generalized to uq and , respectively, as defined in eq. (p2D. 

We now turn to the problem of subtracting the contributions from excited and continuum 
states in sum rules. In the usual approximation based on duality, one identifies the spectral 
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functions p*^'^-^ and p'* beyond a given boundary in the (si, S2)-plane. Then, the subtraction 
effectively amounts to restricting the dispersion integrals in eq. ( PT| ) to the region below this 
boundary. Ideally, the result should not depend on the precise shape of this region. To be 
specific, one may take 

+ 4 < (BIO) 

where sq plays the role of an effective threshold. Popular choices of the duality region are 
triangles in the (si, S2)— plane corresponding to a = 1, and squares corresponding to a — > 00. 
Since the spectral density (P5| ) vanishes everywhere except at si = S2, it is actually irrelevant 
which form of the boundary we adopt i. 

Using duality as outlined above we have to evaluate the integral in eq. ( pT]) with the 
integration region restricted by eq. ( P10|) . Changing variables and integrating over v one 
obtains an expression similar to eq. ([B8|) , but with the upper limit of integration in s lowered 
to 2so and with the addition of surface terms. The latter disappear for = M|. Hence, one 
is again led to eq. (PP|) with a simple modification of the integration limit: 



F(M2) = mj^ 



E 

k 



2k+^ 



2so 



2m2 



— e aJ2 



(Bll) 



This proves the substitution rule quoted in eq. (^). It is important to note that the propor- 
tionality of the Borel transform F{Mf, M|) given in eq. ( P8| ) to the wave function yj^ at the 



point Mo = Mil {Ml + M|) is generally destroyed by the continuum subtraction. It is only 
retained in the symmetric point Ml = M| implying uq = 1/2. 

The above procedure is not possible for higher twist contributions which are proportional 
to zero or negative powers of the Borel parameters. The reason is that the corresponding 
spectral densities are not concentrated near the diagonal Si = S2- In fact, the continuum 
subtraction is rather complicated in these cases. For further discussion we refer the reader 
to the second paper of ref . . Here, we neglect the continuum subtraction in higher twist 
terms altogether. This is justified to a good approximation since the corresponding spectral 
densities decrease fast with si and S2 as a consequence of ultraviolet convergence and, hence, 
the continuum contribution is expected to be small anyway. 



^This is literally true only if the power series defining the wave function (B4) is truncated at some finite 
order, or if it converges rapidly. However, this condition is always met at a sufficiently high normalization 
point where the wave function deviates little from the asymptotic form. 
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Table 1: Summary of theoretical estimates. 



Reference g 93* Bn go'Dw T{D*+ D^TT+) (keV) 

This paper 0.32 ± 0.02 29 ± 3 12.5 ± 1.0 32 ± 5 

This paper'^ - 28 ± 6 11 ± 2 

[|n|]" - 32 ± 6 - 

0" 0.39 ± 0.16 20 ± 4 9 ± 1 

[0]"* 0.21 ± 0.06 15 ±4 7±1 10 ±3 

0.7 - - 

- 64 - - 

0.75 ^ 1.0 - - 100 ~ 180 

3|]^ 0.6 ^ 0.7 - - 61 ^ 78 

0.4 ^ 0.7 - 



0.3 



- - 16.2 53.4 

ra-^ - - 19.5 ± 1.0 76 ± 7 



MP - - 8.9 16 



[H]'^ - - 8.2 13.8 

Experiment' - - < 21 < 89 

" QCD sum rules in external axial field or soft pion limit. 
* including perturbative correction to the heavy meson decay constants. 
Quark model + chiral HQET. 

Chiral HQET with experimental constraints on D* decays. 

Extended NJL model + chiral HQET . 

Quark Model + scaling relation. 
^ Relativistic quark model. 
^ SU(4) symmetry. 
^ Reggeon quark-gluon string model. 

^ Combination of ACCMOR [33] and CLEO [34] measurements. 
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THE FIGURES 





Figure 1: Diagrams contributing to the correlation function (3). Solid lines represent quarks, 
dashed lines gluons, wavy lines are external currents, and the ovals denote tt meson wave 
functions on the light-cone. 



28 




2 2.5 3 3.5 4 




5 6 7 8 9 10 11 12 13 



Figure 2: The r.h.s. of the sum rule (44) for the couphng constants (a) gD'Ov and (b) §3*3^ 
as a function of the Borel mass squared. The arrows indicate the interval in Af^ allowed by 
the reliability criteria specified in the text. 
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Figure 3: The form factors for the transitions (a) D — )• tt and (b) i? — )• tt as predicted by the 
hght-cone sum rule (sohd hues) in comparison to the single-pole approximation (dashed lines) 
with the normalization fixed by the coupling constants gD'Ov and gB*BTn respectively. 
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